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Abstract 

We study short-range ferromagnetic models residing on planar manifolds with 
global negative curvature. We show that the local metric properties of the em- 
bedding surface induce droplet formation from the boundary, resulting in the 
stability of a Griffiths phase at a temperature lower than that of the bulk tran- 
sition. We propose that this behavior is independent of order parameter and 
hyper lattice specifics, and thus is universal for such non-Euclidean spin mod- 
els. Their temperature-curvature phase diagrams are characterized by two 
distinct bulk and boundary transitions; each has mean-field critical behav- 
ior and a finite correlation length related to the curvature of the embedding 
surface. The implications for experiments on superconducting hyperlattice 
networks are also discussed. 
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I. INTRODUCTION 



Disordered spin systems behave in ways quahtatively distinct from their periodic 
counterparts.llhi More specifically, they display broad relaxation spectra, in clear contrast 
to the Debye relaxation observed for spin crystalline materials. Relaxation times depend on 
local conditions, and thus it is not surprising that randomly-coupled spins relax on a wide 
distribution of time-scales. There are now several anomalous magnetic materialJi^ whose 
non-Debye relaxation seems to be determined by lattice topologyB rather than by disorder. 
The theoretical challenge is to identify and characterize periodic reductionist models that 
exhibit such behavior in the absence of intrinsic randomness. 

There are several analytic studies of periodic glass models,! though their infinite-range 
nature makes their relevance to real materials unclear. More specifically, the absence of any 
length-scale in these approaches implies relaxation on just a few time-scales. Detailed ana- 
lytic methods to treat short-range glasses, even with intrinsic disorder, remain to be found. 
It has been proposed that their slow dynamics result from rare, locally-ordered spatial re- 
gions that are energetically probable due to random spin-spin couplings,&0 similar to the 
situation associated with Griffiths phaseJi^ in diluted ferromagnets; here large fluctuating 
droplets lead to relaxation on long time-scales. Though the application of this droplet sce- 
nario to experimental spin glasses remains controversial,lii it is an appealing starting point 
for geometrically-induced glassiness. In quenched ferromagnets, the random initial configu- 
ration approaches equilibrium via domain wall motion.El This coarsening occurs slower when 
the boundaries are pinned by impurities, as is the case in the random-field Ising modehlii 
Thus one can ask whether local geometry can induce such slow domain wall dynamics for a 
short-range ferromagnetic model. Self-similar lattices are excellent candidates, since all mi- 
nority droplets will be non-compact indeed non-trivial slow relaxation has been reported 
in this case.lli'0 It is therefore natural to continue this program by studying ferromagnets on 
lattices embedded in surfaces with constant negative curvature. Dynamics on such hyper- 
bolic manifolds, particularly in the area of chaos,^'^! are known to be qualitatively similar 
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to those observed in disordered system^2l and thus such surfaces provide promising settings 
for the identification of broad relaxation spectra in disorder-free models. 

In this paper we show that periodic ferromagnetic models on hyperlattices with open 
boundary conditions exhibit slow dynamics on a distribution of time-scales in the absence 
of intrinsic randomness. Their metric structure makes it energetically probable for large 
domains of minority spins to nucleate at the boundaries, leading to the stability of a Grif- 
fiths phase at a temperature lower than the bulk temperature. In dilute ferromagnets 
with correlated disorder, the presence of such rare droplets leads to a diverging magnetic 
susceptibility,!! a phenomenon that we also observe for the periodic hyperlattice models. 
We attribute this behavior to the distinction between bulk and boundary sites, where the 
latter comprise a significant proportion of the total site number. Two mean-field transitions 
are identified, associated with the bulk and the boundary respectively, whose characters are 
determined by local metric properties; for example the correlation length associated with 
each transition is related to the curvature of the embedding surface. We conjecture that 
this phase behavior is universal for all short-range ferromagnetic models residing on lat- 
tices embedded on hyperbolic planes, independent of details associated with the spin order 
parameter or lattice specifications. 

The outline of this paper is as follows. In the next Section (Section II) we discuss 
geometrical properties of lattices embedded in surfaces of negative curvature. The model 
and its specific heat are presented in Section III. The Bethe-Peierls transition for the central 
spin, that which is deepest in the hyperlattice, is derived in Section IV. This analysis is 
generalized in Section V to an arbitrary lattice site. Analysis of magnetization distributions 
is presented in Section VI, with particular emphasis on the inferred slow dynamics. In 
Section VII, we summarize our results with a conjectured phase diagram for these models 
in the curvature-temperature plane and discuss possible experiments on Josephson junction 
arrays. In order to maintain the fiow of the text, we have delegated the derivation of the 
main analytical results to Appendices. 
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II. HYPERLATTICES: REGULAR LATTICES ON A HYPERBOLOID 



Infinite regular lattices are characterized by two integers {p, q) where p is the number of 
polygon edges and q is the number of polygons around a given vertex. On an Euclidean plane 
(i.e. no curvature), there are three tilings corresponding to {p — 2){q — 2) = 4: the square 
(4, 4), the triangular (3, 6), and the the honeycomb (6, 3) lattices. On the sphere ^2 (positive 
curvature) the five platonic polyhedra correspond to the condition {p — 2)(g — 2) < 4: the 
tetrahedron (3,3), the cube (4,3), the octahedron (3,4), the dodecahedron (5,3), and the 
icosahedron (3,5). Hyperlattices, tilings of the hyperbolic plane with negative curvature, 
correspond to the condition {p — 2){q — 2) > 4; clearly there exist an infinite number of 
possibilities. The site centered (3, 7) hyperlattice is shown on Fig. |l[ 

Loopless trees, such as Bethe lattices, correspond to regular tilings of the hyperbolic 

plane labeled by the integers (cxD,g).0il Hyperlattices can be built in a shell structure, 

where each layer is analogous to a generation of a tree. For a hyperlattice of n shells, the 

i-i 

number of boundary sites, 6Vn, scales with the number of bulk sites, Vn, since 6Vn ^ Vn 
for a Euclidean tiling, the hyperlattices are considered to have d = oo. However they possess 
intrinsic length scales determined by their radii of curvature, a feature that influences the 
critical behavior of models residing on these surfaces.il 

It is useful to study the continuous-space analogue of the hyperlattice, represented as a 
unit disc on the Lobachevskii plane with a metric 

9{z) = (1) 

(1 - |zP) 

where z = x + iy is a point on the disc. A finite-size system is obtained by restricting z to 
the disc |2;| < -R < 1. The minimum- length paths, e.g. the geodesies, between two points at 
coordinates Zi and Z2 are circles perpendicular to the unit circle. The distance between zi 
and Z2 is 



d{zi, Z2) = tanh 



zi - Z2 
1 - 2:1^2 



(2) 
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FIGURES 



FIG. 1. The site centered (3-7) hyperlattice with 3 generations. 

In order to illustrate the phenomenology associated with different curvatures, we deter- 
mine the short-scale corrections to the perimeter of a circle, £, of radius d for the negative 
and positive cases. First we consider the circle embedded on a two-dimensional sphere where 
its perimeter is decreased by the presence of a positive curvature TZ: 

Cs^ = 27rd-^d^ + .... (3) 

We can perform an analogous expansion for the case of the hyperbolic plane with the metric 
Eq. |I[ The set of points at a fixed distance d away from a given point X (chosen on the real 
axis) is found to be represented by the circle z = xq + Rq exp (i0), with 

1-tanh^ci 2n tanhd 
Xq = X and Ro = (1- X^) (4) 

The periphery of the circle z = Xq + Rq exp (i0) is 

-^^buik = ^ \Jg{z)RQd(t) = TT sinh {2d). (5) 

so that a short-scale expansion to order d^ leads to 
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4lTi 

Cbni^ = 2nd + —d^ + ... (6) 

Comparing Eq. ^ with Eq. ^, we see that the curvature associated with the metric Eq. ^ is 
TZ = —1/2. On large length-scales, the sphere is compact; there are no more points at a 
distance larger than a critical value. By contrast, on the manifold with negative curvature, 
the number of points at a distance d away from a given point increases exponentially with 
d above the curvature radius: 

TT 

>Cbuik ~ -exp (2o?), (7) 

which is obtained from Eq. |^. On a Bethe lattice, the number of points at a given distance 
d scales like z'^, with z the branching ratio. This indicates a close link between the tree and 
hyperlattice structures, related to the underlying manifold with negative curvature. 

III. THE FREE ENERGY OF THE HYPERLATTICE MODEL 

Here we study the nearest-neighbor ferromagnetic Ising model (FIM) with Hamiltonian 
Ti = <^i<^j, where a = ±1 and are neighboring sites on a hyperlattice. A char- 

acterization of a special case, the FIM on the Cayley tree, has already been reported. 
However there exists a temperature scale, Tg, below which large droplets of flipped spins 
proliferate from the boundaries, resulting in non-Gaussian magnetization and glassy behav- 
ior for a macroscopic number of sites. In this paper we investigate such boundary-induced 
Griffiths phases for periodic spin models on more general hyperlattices. 

The free energy per spin of the ferromagnetic Ising model on the Cayley tree is analytic for 
all temperatures; because of the absence of loops on this pseudo-lattice, it can be obtained 
from a high-temperature series expansion and is = — 2 In [cosh(/3)]. The situation 

could be different for a general hyperlattice due to the presence of loops. The number of 
independent cycles tIc of a graph with sites and of bonds is ric = nt — {us — 1). For 
the (3,7) hyperlattice, we have ric/ns ~ (5 — a/5)/2 1.38197, which is even larger than 
the square lattice value ric/ns =~. Thus the free energy of the FIM on the hyperlattice 



could have a singularity, and its absence/presence must be checked explicitly. In order to do 
this, we have measured its internal energy (e) and specific heat = ((e^) — (e)^)/T^ using 
numerical simulations. We have ensured that the equilibrium distribution was properly 
sampled by comparing the numerical estimate of d{e)/dT and ((e^) — (e)^)/T^ and have 
used both heat-bath and cluster algorithms^, with perfect agreement. The temperature- 
dependence of the specific heat is shown on Fig. ^ for an increasing number of layers n = 
8, 11, 12. All these curves display a maximum around the same temperature, but these peaks 
do not increase with the system size; thus we interpret this maximum as a Schottky anomaly 
rather as a signature of a phase transition in the thermodynamic limit. The free energy of 
the FIM on the hyperlattice thus appears to be analytic at all temperatures, similar to the 
situation for the special case of the Cayley tree. We note that here we are always considering 
the case of open boundaries. 
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FIG. 2. Specific heat for a (3,7) hyperlattice computed from the numerical derivative of the 
internal energy and from the fluctuation of energy. The agreement indicates that thermal equi- 
Hbrium is properly sampled. The results presented are obtained with 8, 11 and 12 layers. The 
number of sites is 4264 (8 generations), 76616 (11 generations), and 200593 (12 generations). The 
number of bonds is 10150 (8 generations), 182490 (11 generations), and 477799 (12 generations). 
No sign of divergence is found when the system size is increased. 
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IV. THE BETHE-PEIERLS TRANSITION OF THE CENTRAL SPIN 



In order to characterize the behavior of the FIM on a hyperlattice, we begin by consid- 
ering the ordering of its central spin, that residing on the site deepest in the lattice (e.g. 
the site Tq on Fig. §). As a point of reference, we review the special case of the Cayley tree 
where the central spin is known to undergo a mean-field Bethe-Peierls transition.iH^. In 
response to an applied uniform field on the entire tree, its central spin develops the local 
susceptibility Xo'^'^ 

= E exp = E - 1)" exp . (8) 

Here we use the subscript "0" in Xo'^'^ to emphasize that we are only considering the be- 
havior of the central spin. In Eq. |^, di is the distance between the central site Tq and the 
site /; = — 1/ In (tanh (/?J)) is the correlation length set by the exponential decay of 
the spin correlations {<Jk<Ji) ~ exp {—dk,i/^T), identical to that of the Ising chain because 
of the absence of loops. The central spin susceptibility x^o'^'^ diverges if T < Tbp, with 
{z — 1) tanh (/^Bpt/) = 1. We emphasize that this transition results from the metric of the 
embedding curved space, which leads to the to the prefactor {z — 1)" in Eq. this is the 
number of sites at generation n from the central spin. Thus, thanks to this prefactor, expo- 
nentially decaying correlations are sufficient to generate mean-field ordering of the central 
spin. As an aside, we note that here we have only considered x^o'^'^^ but this calculation can 
be generalized to characterize the full thermodynamic behavior associated with this spin.!!!! 

For the FIM on the general hyperlattice with loops, the expression for the local suscep- 
tibility is also dominated by the long length-scale contribution of the spin-spin correlations, 
leading to a diverging susceptibility with a finite correlation length.ii Eq. |^ may be adapted 
to evaluate the central-spin susceptibility using a continuous-space description 



Tx'"\0)= f g{z)exp(-^]d'z, 
J\z\<R \ ) 



(9) 



where g{z) is the metric in Eq. |I|, and d{z^ 0) is the distance between the origin and the 
site aX z = pexp{i9) (see Eq. ||). We implicitly identified the distance on the graph and 
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the "hyperbolic" distance given by the metrics. In fact, for a given distance on the graph, 
there is a Gaussian distribution of hyperbohc distances, which has been studied in details 
for a Cayley tree model by Comtet, Nechaev, and Voituriez@. We claim that the physics of 
the FIM is not sensitive to the existing difference between the two distances. The resulting 
susceptibility 



diverges when C,t = 1/2? with the mean-field exponent 7 = 1. The correlation-length 
exponent can be obtained by matching the Euclidean and hyperbolic spin correlations on 
short length-scales at the Bethe-Peierls transition. The Euclidean behavior is (cxcry) ~ 
l/|x — y|'^^^+'' with d = 2, while the spin correlations decay exponentially on the hyperbolic 
plane; this argument leads to r] = 0. 

We have therefore identified two critical exponents, 7 and r], to have mean-field values. 
Scaling relations imply that all exponents of the hyperlattice Bethe-Peierls transition are 
then mean-field as well. We note that there is another way to obtain this result as well. On 
(i-dimensional Euclidean lattices. The upper critical dimension of spin models is related to 
the probability of intersecting random walks.0 the upper critical dimension of spin models 
on a (i-dimensional Euclidean lattice is related to intersection properties of random walksS. 
On hyperlattices (including the special case of trees), the return probability of a random 
walk is vanishingly small@H^, indicating that FIM hyperlattice models are above their upper 
critical dimension, and thus mean-field transitions are expected. 
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FIG. 3. A tree with 4 generations and a forward branching 2 — 1 = 2. We denote by Tq ••• 
the vertices along a path from the top site to a leave site. The sites descendant from a given 
vertex 7^ are grouped together as shown on the figure to calculate the generalized Bethe-Peierls 
susceptibility in section V. 



V. GENERALIZED BETHE-PEIERLS TRANSITIONS 

We would like to consider a possible Bethe-Peierls transition for an arbitrary spin of the 
FIM on a general hyperlattice. Again, we begin by studying the special case of loopless 
Cayley trees. Here we consider spin ordering on an arbitrary site 7^ at generation m = 0, 
...,77. (see Fig. ^. We group the sites, indicated in Fig. ^, to expand the susceptibility in 
powers of a; = tanh (/3J) (for details see Appendix ^ to obtain 

Tv*- _ 1 - [(^ - 1)0:]—+^ 1^ 

^™ " l-(^-l)x ^ \-x ^ ^ 



- 2)x2 
^ 1 - iz-Dx 



— X 1 — (2; — \)X^ 

In the limit n,m ^ +00 and n — m constant (boundary behavior), x^m'^{T) is smooth at 
the Bethe-Peierls transition temperature {{z — 1) tanh (PbpJ) = 1), but diverges at a lower 
temperature T' corresponding to (-2 — 1) tanh^ iP'J) = 1- As shown explicitly in Appendix 0, 
the existence of this transition results from a combination of the local metric (i.e. the 
number of sites a distance from the reference one) and exponential spin correlations. In the 
limit n —* +00 and m fixed (bulk behavior), x^rn"^ i^i Eq. |10| diverges at the Bethe-Peierls 
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temperature Tbp with mean-field behavior. If n +00 and m = An, the susceptibility 
in Eq. |1^ shows a divergence at the Bethe-Peierls temperature. The temperature T' was 
already found by several authors by studying the magnetic field dependence of the Potts 
model on the Cayley tree@0. Here, we provide another interpretation for this transition, in 
terms of a local Bethe-Peierls transition controlled by the metric properties at the boundary, 
and prove that this transition is related to a Griffiths phase. 

We have therefore demonstrated the existence two distinct Bethe-Peierls transitions as- 
sociated with the bulk (where n — m +00), and the boundary, related to the different 
metric properties associated with these two different site species. They also affect percolation 
thresholds on these trees, as discussed in Appendix 0. 

Now that we have understood the bulk and boundary transitions of the tree FIM in 
terms of metric properties, we continue to ask a similar question about their more general 
hyperlattice analogues. We stress that on these manifolds of negative curvature, site-scaling 
is different for reference sites in the bulk and on the boundary. As displayed in Eq. |^, 
the number of sites at a distance d from a given bulk reference one scales like £buik ~ 
(7r/2) exp {2d) in the continuous space model. In Appendix ^ we show the analogous quantity 
associated with a boundary reference site scales like C ~ exp (d). This result is compatible 
with our numerical determination of the site-scaling at a given distance, as displayed in 
Fig. |. 
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(3,7) hyperlattice with 14 generations 




5 10 15 20 25 30 

k 

FIG. 4. Logarithm of the number of sites of the (3, 7) hyperlattice with fourteen generations as 
a function of distance k from a given site located (i) on the center of the hyperlattice (diamonds) 
where the calculated exponential growth is indicated by a dashed line (ii) on the boundary (crosses) 
with calculated exponential growth (short dashed line) that is the square-root of the calculated 
bulk growth (dashed line) (iii) on an intermediate (7) generation (boxes), where a crossover is 
observed between bulk and boundary behavior with increasing distance scales. 

The expression for the local susceptibility involves a correlation length, which is not 
necessarily uniform for the full lattice. We assume a uniform, isotropic correlation length as 
an Ansatz and then deduce the existence of bulk and boundary transitions for the general 
hyperlattice set by = 1/2 (see section [TVl ) and by = 1- The existence of two different 
transitions does not rely on the initial Ansatz because the correlation length at the boundary 
can only be reduced compared to that of the bulk. Therefore, the initial Ansatz of uniform 
correlations indicates the existence of a boundary transition at a temperature lower than 
the bulk transition. 

We support the previous analysis, assuming a uniform spin correlation length, with 
numerical determinations of the spin-spin correlations. On the Cayley tree, the correlation 
of two spins a and a' at a distance d is [tanh {pd)Y, irrespective of the site locations of 
the two reference spins. By contrast, as discussed below, the situation is different for the 
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general hyperlattice case where the correlations are reduced at the boundary. We distinguish 
between spin correlations in the bulk and at the boundary: 

Bulk: The symmetry group of the (3, 7) hyperlattice with a finite number of layers n 
is generated by one rotation of angle 2^/7 around the central spin, and one inversion with 
respect to an axis going through the central spin. The infinite hyperlattice has a huge 
symmetry group: the lattice is left invariant under the aforementioned set of symmetries 
around any lattice site, a fact that is no longer true for a finite system. However the 
properties of the finite-size lattice refiect the huge symmetry of the infinite lattice. To be 
more precise, let us consider a FIM on a hyperlattice with n layers and Nk correlations 
{<^o<^k) between a given spin (Tq, and any of the A^^ spins at a distance k from ctq. Let us first 
take (To to be the central spin. Strictly speaking, the number of such different correlations 
is Nk/7 because of the 27r/7 rotation symmetry. However we find strong evidence that, 
for k fixed and n large, all the (aocrk) correlations nearly coincide. On the lower, negative 
X-axis of Fig. ^, we have chosen k = 2 and increased n; the effect of the site alternation 
with coordination 3 and 4 is clearly visible (see the boundary on Fig. |I|). The correlation 
inhomogeneities are reduced with increased system size (cf. the negative x-axis part of 
Figs. |and§). 

Boundary behavior: Now we consider the reference spin ctq to be located at the hyper- 
lattice boundary and the spins cr^ at a distance k from ctq. The dispersion is similar to that 
already described for bulk case, though it persists in the large-n limit (see the positive x-axis 
part on Fig. ^). The sites with smaller coordination are less correlated, as expected on phys- 
ical grounds. Therefore the correlations are weakened at the boundary of the hyperlattice. 
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FIG. 5. Spin-spin correlations at T = 1 for the central (negative x-axis) and a boundary spin 

(positive X-axis) at a distance k = 2 where x is an arbitrary spin index. 



o 

I 

o 
o 



(3,7) hyperlattice with n=ll generations 
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FIG. 6. Spin-spin correlations computed for the central and a boundary spin at a distance 



k = 7. 



VI. MAGNETIZATION DISTRIBUTIONS AT LOW TEMPERATURES 

We now demonstrate slow dynamics at low temperatures in this periodic hyperlattice 
FIM. Towards this goal,we track the energy and magnetization distributions, comparing 
results obtained using two sampling methods. The first is a standard heat bath single 
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spin-flip algorithm that probes configuration space where two configurations are considered 
neighbors if they differ by exactly one spin. The second one is a cluster algorithm that 
can go from any configuration to any other in just one step, and is applicable here because 
of the absence of frustration. We use these two algorithms to probe the structure of the 
configuration space, analogous to similar studies performed previously on Cayley trees.liiiil 
In both cases, the valleys correspond to magnetic domains in real space (see Fig. |^). On 
the hyperlattice, the energy of a magnetic domain scales like the logarithm of its area (see 
Appendix 0) while it does not scale with area on the tree. 




(a) (b) 

FIG. 7. Magnetic domains (a) on the tree, and (b) on the hyperlattice. The geodesic 1-3-2 
is perpendicular to the border of the hyperlattice Rexp Point 1 is Rexp{i9). Point 2 is 

Rexp (—iO) 

In Figure ^, we present three energy distributions. One was acquired with a cluster 
algorithm, while the other two were obtained with the single-spin fiip algorithm and 3 x 10^ 
and 4.78 x 10^ MCS-spins. The three simulations were performed at a temperature T = 1.4, 
below the estimated Bethe-Peierls transition (of order of the lattice coordination). These 
three energy distributions are close to the same Gaussian distribution, centered at the mean 
value of the internal energy, with a width proportional to the specific heat times the square 
of the temperature. iFiom this figure it is clear that there are no long time-scales present 
as far as the energy is concerned. 
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However, the behavior of the magnetization is rather different. In Fig. |] we present 
three magnetization distributions, determined during the same runs as the internal energy 
distributions. The three magnetization distributions are distinctive. The two distributions 
obtained using the single-spin flip algorithm are bimodal whereas the distribution obtained 
with a cluster algorithm is unimodal. For other runs, with less statistics, the magnetization 
distribution with the single spin flip algorithm was found to be multi-modal, with three 
or more maxima. This suggests the existence of energy barriers which trap the system 
with single-spin dynamics. These simulations provide evidence that the low temperature 
dynamics are slow. The glass crossover temperature corresponds to the correlation length 
being comparable to the system size. Also the magnetization distributions are broad at low 
temperature (see Fig. ^ for the hyperlattice, similar to that found for the Cayley tree;0 
this is an indication of relaxation on many time-scales. This is because of the stability of 
a Griffiths phase on the entire lattice, at a temperature set by the boundary Bethe-Peierls 
transition (see Appendix ^). 

Energy distribution at T=1 .4 for a (3-7) liyperlattice with 1 1 layers 
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FIG. 8. Energy distribution computed at T = 1.4 using single-spin flip cluster algorithms where 
"more" and "less" curves correspond to 47800 and 300000 MCS/spins respectively; the "cluster" 
curve was obtained with 10200 updates of the cluster algorithm. 
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Magnetization distribution at T=1 .4 for a (3-7) hyperlattice with 1 1 iayers 
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FIG. 9. Magnetization distribution computed at T = 1.4 using single-spin flip and cluster algo- 
rithms where "more" and "less" curves correspond to 47800 and 300000 MCS/spins respectively; 
the "cluster" curve was obtained with 10200 updates of the cluster algorithm. 



VII. DISCUSSION 

In summary, we have studied the nearest-neighbor ferromagnetic Ising model on tilings 
embedded in hyperbolic surfaces with negative curvature, often beginning with the spe- 
cial case of Cayley trees. We have identified two mean-field transitions in these systems, 
associated with ordering of the bulk and the boundary spins as a function of decreasing tem- 
perature. These two transitions can be understood in terms of the local metric properties 
of the lattice, specifically by the distinct scaling of sites with distance from the boundary 
and bulk spins. We beheve that these two transitions will be characteristic of any short- 
range ferromagnetic model residing on a planar surface with globally negative curvature 
independent of lattice specifics or the nature of the spin order parameter. 

In order to satisfy the criterion of "global negative curvature" , at large length-scales the 
resident lattice should be embedded on a surface with a negative (not necessarily constant) 
curvature instead of on a Euclidean plane. For instance, we could imagine models with a 
fluctuating metric. As a simple example, we could consider a recursive random lattice in 
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which the sites have a coordination z with a given probabihty Pi^z). The sites with z = 2 
have the local environment of an Ising chain, and the sites with 2; > 3 have local branching. 
The bulk Bethe-Peierls transition is found to be ((2; — 1)) tanh (/5buik</) = 1, whereas the 



boundary transition corresponds to ^((^ — 1)) tanh (/^bound^) = 1- Both transitions occur 
at a finite temperature as soon as ((2;)) > 2, in which case the metric is globally hyperbolic. 

One can also question how the hyperbolic/Euclidean transition would operate in a con- 
tinuous model in the limit of a large curvature radius 7?. ^ 1. The bulk and boundary 
Bethe-Peierls transitions would be of order 1/7?. while Euclidean physics would develop be- 
low the length scale 7?.. Strictly speaking, the Euclidean phase transition occurs only when 
Tt = +00. However, with a finite but large TZ, we expect a pronounced maximum in the 



susceptibility and specific heat. The expected behavior is summarized on Fig. |Ty. Again 
we emphasize that we expect this phase behavior for ferromagnetic models residing on hy- 
perboloids independent of the details of the spin order parameter. For Cayley trees, it has 
been shown that the nearest-neighbor xy model displays similar physics (i.e. mean-field 
transition) as its Ising counterpart.0 For a general hyperlattice, we expect two mean-field 
transitions for short-range xy models since there exist loops and two species of spin (bound- 
ary and bulk), as in the Ising case. As was done in this Paper, these transitions would 
be obtained by comparing the exponential decay of the correlations to the growth in the 
number of sites with distance; this would be interesting to verify.. 
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FIG. 10. (Curvature, Temperature) phase diagram. Tc is the critical temperature of the Eu- 
chdean model. The dashed line is a cross-over. Region (1) is the paramagnet. Region (11) corre- 
sponds to algebraic correlations below the curvature radius TZ, and exponential correlations above 
TZ. The solid lines are the bulk and boundary Bethe-Peierls transitions, determined as = ^ 
and = 27^ respectively. Note that (^t is related to the exponential decay of the correlations 
above the curvature radius, and not to the behavior of the correlations below the curvature radius. 
Region (IV) corresponds to the Griffiths phase of the entire lattice. 

The physics of the low-temperature behavior of the hyperlattice FIM is determined by the 
formation of droplet-hke excitations nucleating from the boundary. These rare fluctuations 
do not affect the specific heat (Section 111), but do lead to a broad magnetization distribu- 
tion reflecting a wide spectrum of relaxation time-scales. Such Griffiths phases are usually 
associated with dilute ferromagnets with spatial inhomogeneity. In the periodic models we 
have studied here there is no intrinsic disorder; however there are two distinct site species, 
a feature that contributes to favoring droplet formation. Indeed here we report a diverging 
susceptibility per site at the boundary transition(see Appendix |^), a characteristic shared 
with random ferromagnets with correlated disorder^. Similarly divergences are reported for 
quantum disordered magnets, where the randomness is correlated in the time-dimension. 
To our knowledge, this is the flrst time a Griffiths phase have been identifled in a periodic 
magnetic system with loops; furthermore the identiflcation of slow dynamics in a short-range 
periodic system at low temperatures is very encouraging. 

Though motivated by experiment, our chosen lattices of study are somewhat removed 
from the structures commonly observed in Nature. However perhaps we can close by at- 
tempting to reconnect with experiment. In particular, we have identifled Griffiths phases 
in ferromagnetic models residing on hyperlattices. Though this phenomenon has been dis- 
cussed extensively in the theoretical literature, it has not yet been conclusively identifled 
in the laboratory.^! For example, there have been claims in random-fleld materials,^ but 
they remain controversial due to plausible alternative interpretations of the data.0 Perhaps 
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Josephson junction arrays fabricated in a hyperlattice topology would provide a promising 
setting for the observation of this well-discussed phenomenon. In this artificial network, a 
junction is located at each link of the hyperlattice; the short superconducting wires have a 
phase (pi with a Hamiltonian 

H = EjY,[1-cos (0.-0,)]. 

where the ground-state configurations corresponds to a uniform 0j = and Ej is an energy- 
scale associated with a junction. If this array were placed in a time- dependent transverse 
magnetic field, the low-frequency part of the resulting ac susceptibilityil should provide a 
useful probe for the slow dynamics associated with the Griffiths phase. 

In order for our ideas to be applicable to the Josephson hyperlattice, it should not 
display a vortex binding-unbinding transition. Assuming that the nature of the spin order 
parameter is not important, we can apply some of our previous expressions to check this. 
The phase at a distance d away from a given bulk vortex scales like 0buik(O ~ 4exp (—2/), 
where we have used the expression Eq. for the number of sites at a given distance I. 
The energy of a junction a distance d away from a bulk node leading to the 

total vortex energy -Ebuik = '^t^Ej. It is remarkable that this quantity is finite, since it 
diverges logarithmically with distance on the square lattice. Similarly, the energy of a 
vortex at the boundary is found to be -Ebound = 27r^i?j, again a finite quantity. As a 
consequence, on the hyperlattice, there is no Kosterlitz-Thouless transition associated with 
the unbinding of vortex-antivortex pairs. We note that an identical conclusion was reached 
for two-dimensional Coulomb systems residing on a surface of constant negative curvature;^ 
this agreement gives us further confidence in our conjecture about the universal nature of 
the phase behavior that we have identified. 
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APPENDIX A: NUMBER OF SITES AT A GIVEN DISTANCE ON THE TREE 



We consider a site at generation m on a tree and calculate the number of sites Nm{l) at 
a given distance / from the site at generation m. This is given by the prefactor of the term 
in the expansion of the susceptibility in powers of x = tanh (/3J): 

Tx'^pira) = ^ - - ^)^r-"'^' (Ai) 



l-{z- l)x 

+ x"^'^ fl + (-2 - 2)x + {z- 2){z - l)x^ + ... + {z- 2){z - 1) 



fc=0 

Eq. |iy has been obtained by summing the series in Eq. First the sum Y^1Z^[{z — 
gives rise to {z — 1)' sites at a distance < / < n — m. Next, the sum Y1^=q x"^~'' results 
in one site at a distance 1 < I < m. The remaining terms give rise to {z — 2)(z — 
sites at a distance /, with < k < m — l and m — k + l<l<n + m — 2k. We therefore need 
to distinguish between two cases: case (I) m > [(n — l)/2]; and case (II) m < [{n — l)/2], 
with [...] the integer part. 

• Case (I): m > [(?7, — l)/2]. We should compare I with n — m and m + 1. We find three 
cases: 

— {i) I < n — m. The number of sites at distance I is {z — 1)' — 1. 

— (ii) n — m<l<m + l. The number of sites at distance I is [z — l)[("~'"+')/2] _ i. 

— (iii) m + 1 < /. The number of sites at distance I is (2- 1)[("-™+0/2] _ 

• Case (II): m < [(n — l)/2]. We compare again / with n — m and m+1: 

— (i) Z < m + 1. The number of sites at distance / is (z — 1)' — 1. 

— (ii) m + l</<'n, — m. The number of sites at distance / is (2; — 1)' — {z — 1)'^™"^^. 

— (iii) n — m < I. The number of sites at distance / is (2; — 1)[("~™'+')/^] — (2; — 1)'^™"-*^. 

Therefore at large distance the number of sites at a distance / from a given site at generation 
m scales like A^m(0 ^ {z — 1)'''^ if m > [{n — l)/2]. This scaling should be compared with 
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the exponential decay of the correlations {<7i<7j) ~ exp (— (iij/^r) and leads to a boundary 
transition at a temperature T' given by {z — 1) tanh^ iP'J) = 1- This transition occurs 
whenever n — m is finite in the limit n ^ +00. 



APPENDIX B: PERCOLATION ON THE TREE: GENERALIZED 

BETHE-PEIERLS LIMIT 

We would like to illustrate the different behavior of the generalized Bethe-Peierls tran- 
sitions on the example on bond percolation where the same phenomenon occurs as in the 
case of ferromagnetism. We denote by P^^{M) the probability to find M sites (a mass M) 
in the cluster containing the top site of a tree with n generations. We have 

6*1 6»2_i Ml M^_i i=l i=l \ i=l J 

with 6i = 1 with a probability p{l) = fi and 6i = with a probability p{0) = 1 — fi. The 
average number of sites in the cluster containing the top spin of a tree with n generations 
is iterated as Mm+i = {z — + 1. The fixed point value is M* = — {z — l)/i) and 

diverges at the bulk percolation threshold /ip = 1/(2; — 1). 

We now consider the number of sites M{%n) in a cluster containing the point %n in the 
presence of n generations (see Fig. |^). We find 

^^^-^ - ^^^Tjj;^^ ^^^^ 

{z - 2)/i^ / _ 

of the same form as Eq. |10[ If n, m — +00 with n — m constant, the cluster size containing 
the vertex Tm does not diverge at the bulk percolation threshold /Xp while it diverges at a 
larger percolation threshold /ip = 1/ \/z — 1. If n ^ +00 and n — m growing faster than n, 
we find a percolation transition at the bulk percolation threshold yUp. Therefore, similarly to 
ferromagnetism, percolation shows a different behavior depending on how the Bethe-Peierls 
limit is taken. The boundary percolation threshold /ip is larger than the bulk percolation 
threshold /ip. 
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APPENDIX C: NUMBER OF SITES AT A GIVEN DISTANCE FROM A 
BOUNDARY SITE ON THE HYPERLATTICE 



We consider a finite size disc \z\ < R with the hyperbohc metrics Eq. [l|, and calculate 
the length of the set of points at a distance d from the boundary site at coordinate R. This 
set of points is represented by the circle xq + -Ro^xp (see Eq. which intersects the 
boundary \z\ = R aX the points = xq + iRo exp {±i(po), with 

tanh d 



cos (po 



1 + (-3 + tanh^ d)R^ + (tanh^ d)R^) . 



2R{1 - tanh^ d) 

A straightforward calculation leads to the length of the arc z'^-z~: 



bound 



g{z)Roc 



4 tanh d 
1 — tanh^ d 



■ tan 



-1 



2R 



R^) tanhfi 



\ 2i?+ (1 + i?2)tanhc/ 



In the relevant regime 1 <^ d <^ tanh ^ R with i? ~ 1, we find /Abound ~ exp((i). This 
boundary behavior should be contrasted with the bulk behavior £buik ~ (7i"/2) exp {2d) (see 
Eq.0). 
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APPENDIX D: THE LENGTH OF MAGNETIC DOMAINS 



Unit circle. 





(a) (b) 
FIG. 11. Isometry transform of the droplet excitation of the hyperlattice. (a) the original 

droplet, and (b) the transformed droplet. The isometry is chosen such as the segment l'-3'-2' is on 

the y axis with 3' at the origin. Point 1 is iiexp {iO), point 2 is i?exp {—iO). Point 1' is ip, point 

2' is —ip. 

We show that the length of the droplet excitations shown on Fig. |^(b) is proportional 
to the logarithm of its area. The difference with the calculation in Appendix |C| is that we 
consider here the geodesies 1-3-2 corresponding to the physical situation where the energy 
of the domain wall is minimized. 

To calculate the shaded area on Figs. |^(b) and [ll](a), and use the transformation z — *■ 
z' = {az + b)/ {bz + a) with a and b real numbers, to map the domain on Figs |ll|(a) and ^(b) 
into the domain on Fig. |ll|(b). Imposing point 1 to be transformed into 1' at coordinate ip 
leads to 

^ i?COS^(l+p2) 



a = 



with 



2Rsm9 



a 1 + p'^R^ + 2pRsm9 
1-R^) + [{1-R^f + 4R 



2 ■ 
sm 
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The area element is found to be 

dA = g{r')dr'de' = -7- — ; \ ^ ^ — 7. d^- 

' 2(1 -r2) l + Q;V2 + 2arcos^ 

A straightforward integration over 9 leads to 

' l-aR 



1 + aR 



tan - 
\2. 



When ~ 1, we have A ~ 9/[2{l-R)]. The length of the line 1-3-2 is £ = 



In [(1 + — p)] ~ — I In (1 — i?). The scaling between the area and length of a magnetic 



domain in the limit — > 1 is finally found to be ^ ~ | exp {2C): the area is exponentially 
large in the boundary length. 

APPENDIX E: TOTAL SUSCEPTIBILITY: GRIFFITHS TRANSITION 

We derive the total susceptibility of the tree and hyperlattice models. The suscepti- 
bility of both models diverges below the Griffiths temperature set by the local transition 
temperature of the boundary. 

1. Tree model 

The susceptibihty of the entire lattice Xtot is obtained as the sum of the local suscep- 
tibihties: Xtot = Y.m=oiz — l)™Xm(7'). First if T > T', the susceptibility per spin behaves 
like Xtot /Nn ~ (1 -|-x)^/[l — {z — l)x^], and diverges at the temperature T'. Now if T < T', 
the susceptibility per spin behaves like XTOx/A^n ~ [{z — 1)2;^]"', an infinite quantity in the 
limit n — > +00. 

2. Hyperlattice model 

We first estimate the local susceptibility of the hyperlattice with a finite size \z\ < 
R at the point represented by a real number x, and we assume the correlations to 
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be uniform as discussed in Section V of the main body of the text. We note d = 
d{x,R) = tanh~^ ^ xR)] the distance between the points x and R, and D = 

tanh~^ [{R + a;)/(l + xR)] the distance between the points x and —R. Based on Fig. || and 
Appendix we approximate the growth of the number of sites J\f{l) at a distance / away 
from the point x as follows: (i) I < d (bulk behavior): A/'(/) ^ (7r/2) exp (2/); (ii) d < I < D 
(boundary behavior): A/'(/) ~ (7r/2) exp {I + d). The local susceptibility is found to be 

r-D 

10 



Tx(x) ~ Af{l)exp{-l/^T)dl 
Jo 



TT 1 



2 2 - l/^T 



e fr' — 1 



^ 1 d 

-e 



(El) 



Xtot = / 2'Kxg{x)x{x)dx. (E2) 

JO 



2 1 - l/^T 

In the case of the central spin x = 0, we find a susceptibility diverging at = 1/2. In the 
case of a boundary site x = R, we find a susceptibility diverging at .^t = 1, consistent with 
the behaviors obtained in the main text of this paper. The total susceptibility is obtained 
as the sum of the local susceptibilities Eq. ^: 

rR 

10 

We note R = 1 — c/j, and expand the susceptibility Eq. ^ in the parameter e/j ^ 1. The 
number of sites A^^ scales like e^^. If < 1, the susceptibility per site is found not to 
scale with Ng. If > 1, the susceptibility per site is found to scale like e^^^^''^, an infinite 
quantity in the limit eR 0. 

Therefore, the tree and hyperlattices models show the same behavior: the susceptibility 
per site is infinite below the Griffiths transition temperature T'. The Griffiths temperature 
T' is equal to the temperature of the boundary Bethe-Peierls transition. Physically, this 
originates because of the finite fraction of spins at the boundary: the boundary dominates 
the behavior of the total susceptibility. 
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